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Abstract 

If the universe is slightly non-extensive, and the distribution functions are not exactly given by 
those of Boltzmann-Gibbs, the primordial production of light elements will be non-trivially modified. In 
particular, the principle of detailed balance (PDB), of fundamental importance in the standard analytical 
analysis, is no longer valid, and a non-extensive correction appears. This correction is computed and its 
influence is studied and compared with previous works, where, even when the universe was considered 
as an slightly non-extensive system, the PDB was assumed valid. We analytically track the formation of 
Helium and Deuterium, and study the kind of deviation one could expect from the standard regime. The 
correction to the capture time, the moment in which Deuterium can no longer be substantially photo- 
disintegrated, is also presented. This allows us to take into account the process of the free decay of 
neutrons, which was absent in all previous treatments of the topic. We show that even when considering 
a first (linear) order correction in the quantum distribution functions, the final output on the primordial 
nucleosynthesis yields cannot be reduced to a linear correction in the abundances. We finally obtain 
new bounds upon the non-extensive parameter, both comparing the range of physical viability of the 
theory, and using the latest observational data. 
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1 Introduction 



Primordial nucleosynthesis provides an interesting testing arena where to analyze the viability of physical 
theories, particularly, of the statistical description. It is in this epoch where the earliest bounds upon a given 
theory with cosmological influence can be imposed. Thermal processes (see Ref. Q, hereafter referred as 
Paper I) are non-trivially modified by a non-extensive correction to quantum distribution functions. Then, 
different abundances of light elements are a possible outcome. 

Some of the predictions for primordial nucleosynthesis in a non-extensive setting have been analyzed 
before by some of us, using the asymptotic approach of the quantum distribution functions, see Refs. |@, 
Here, instead, we shall consistently continue within the formalism given in Paper I. Since this approach is 
simpler, we shall obtain analytical results far beyond the point where previous works have reached, see for 
instance Ref. Together with Paper I, we shall then provide a complete history of the early universe, 

accordingly modified taking into account a non-extensive setting. In this paper, we shall focus on the study 
of the changes that non-extensive statistics introduces in the principle of detailed balance, for which we 
provide a detailed analysis, both numerical and analytical. We shall then enhance the study presented in 
1^] , by framing it in a larger picture which encompasses an smaller number of approximations and a larger 
number of predictions. 

Primordial nucleosynthesis was recently used as well to demonstrate that macroscopic samples of neu- 
trinos in thermal equilibrium are indeed distributed according to Fermi-Dirac statistics H]. These latter 
authors considered that neutrinos were distributed by a Bose-Einstein statistics, and established bounds 
(not very restrictive though) upon this unexpected possibility. It is interesting to compare with our case: 
we assume that neutrinos are indeed 1/2 particles, as a large amount of data coming from particles accel- 
erators show, but that even when continuing being fermions, and fulfilling the Pauli's exclusion principle, 
their distribution functions could slightly deviate from an exact Fermi-Dirac one. 

Since we have provided a detailed account of non-extensive statistics, and the reasons why we choose 
the analytical form of the quantum distribution functions we adopted (together with its derivation) in 
Paper I, we shall skip such topics here. We have also considered in Paper I some physical reasons why to 
expect that Boltzmann-Gibbs distribution functions could be considered as an approximation. The same is 
valid for citations to previous works, for which we adopted here the criterion to quote just those needed for 
the explicit computations we are carrying on. This does not mean that our references are the only works 
in cosmological applications of non-extensivity, but only that for proper citation of some of the others, we 
refer the reader to Paper I. 

The layout of this work is as follows. Section 2 presents the basics of the neutron to proton ratio in an 
evolving universe. This section does not use much the fact that we are immersed in a non-extensive setting, 
but just presents general results which are valid disregarding the statistical mechanics used. Indeed, the 
derivation being presented in Section 2 was already given by others |^ |^, and we provide it here just for 
the ease of discussion. In Sections 3 - 7 we give the details of the analysis of the principle of detailed 
balance, and show how to obtain a priori results on the possible range of physically admitted values of 
{q — 1) without the need to compare with experiments. Much of it is done in an analytical form, some is 
solved numerically. In Section 8, we present a detailed comparison between the two situations (full and 
approximate cases) that we found possible for the principle of detailed balance. Towards the end of this 
latter Section we provide a comparison with the latest data available. In Section 9 we compute, for the 
first time in a non-extensive framework, which is the modified capture time, the time in which neutrons are 
captured into deuterons. Using this result we are able to compute the primordial abundance of ^He with 
a greater degree of precision than that obtained in all previous works. We show that there are non-linear 
effects introduced by the appearance of a slight non-extensivity. Finally, we give some general discussion 
in our concluding remarks. 
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2 The neutron to proton ratio 



We begin by turning again to the issue of the evolution of the neutron abundance as the universe evolves. 
We shall base this discussion in the work by, Bernstein, Brown and Feimberg |Q. As we have done before, 
we shall denote by \pn(T{t)) the rate for the weak processes to convert protons into neutrons and by 
Xnp(T{t)) the rate for the reverse ones p. X{T{t)) will be, as usual, the number of neutrons to the total 
number of baryons. For it, a valid kinetic equation is 

^ = A,„(T)(1 - X{t)) - \np{T)X{t). (1) 

The solution to it is given by 

X{T) = f dt'I{t, t')Xpn{t') + X{to)I{t, to). (2) 
J to 

Here, /(t, t') is 

I{t, t') = exp diA{i)^ , (3) 

with 

A{t) = Xpn{t) + Xnp{t). (4) 

Note that this solution is completely general, and does not depend on the statistical mechanics used, except 
by the implicit changes introduced in the new reaction rates. As explained in Q, we simplify by taking 
to = and omitting the X(tQ)I{t,tQ) term. These approximations yield 

X{t) = f dt'l{t,t')Xpn{t'). (5) 
JO 

To compute Eq. (Q), we need to know the reaction rates. Let us consider Xnp{t): 

Xnp = Aj^_|_^^p_|_g- + Ag+„^p_|_j7 + A„^p_|_g-_|_j7 (8) 
that are individually given in Ref.||^: 

roo 

X,+n^p+e- = A / dpAE,{l - (9) 

Jo 

roo 

Xe++n^p+u = AJ^ dpeplp,E,{l - Hf, (10) 

A„^p+g-+p = A r dpepipMi - mi - n, (11) 

Jo 

with A a constant, fixed by the experimental value of A„^p_|_g-_|_p, p^^e are the neutrino and electron 
momenta, and E,^^e their energies. In the energy domain we are interested, some approximations are 
in order, see Refs. [^, ^ for discussion: 1) energy conservation is E^, + rUn = E^ + rup, to be used in 



Finally, we note that 
or, equivalently, 
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Eq. (D and E„ + TUp = Ee + rUn, to be used in Eq.(|oI). 2) In Eq. (0), = Am - Ee > 0, with 
Am = mn — mp = 1.29MeV, from here comes the upper hmit of the integration range. 3) Pauh blocking 
factors are assumed equal to 1. 

Since computations will involve the inverse processes of Eqs. (§) and (|lO|), we shall quote their form 
below: 

A.-^_„^„ = A / dpepipuE.ii - nr, (12) 



POO 

A / dp,pip,E,{i - nr, (13) 

J Am 



where Ei^ + mn = E^ + mp for Eq. (12) and Ei, + mp = E^ + mn for Eq. (|T3|). The lower limit of the 
integral ( |T2|) is given by the minimum momentum that electrons must have in order to result E,y > 0, 
= (Am^ — m^)^/^. On the other hand, the lower limit in the integral (13) comes from the constraint 



E,y = Eg + Am > Am, since we suppose that final states are unbounded. 

As remarked in |^], in general, the electron and neutrino temperatures, Te and T^, may differ because 
at the end of the freezing out period, electrons and positrons annihilate, heating only the photons and 
maintaining with them thermal equilibrium. Indeed, in Paper I we have computed the amount of this 
deviation for our non-extensive setting. This difference is, however, small, and we shall follow Bernstein 
et. al. and set all temperatures equal, T = = = T^. In the standard case this assumption ensure 
that the rates for reverse reactions, such as e~ + p — > n + i^, obey the principle of detailed balance (PDB). 

3 Principle of detailed balance 

The PDB states that if we know the expression for a reaction rate, say Xnp, then this is related with the 
reaction rate for the inverse reaction by an exponential factor, 

Xpn = e-^"^/^A„p. (14) 

We shall study what happen in the non-extensive case, by analyzing each reaction appearing in Eq. (^) in 
a separate way. 

3.1 u + n ^ p + e~ 

We would like to see what is the relationship between A''_ , , and A"^ , , _ , with 

POO 

K-+p^n+. = A , dpepipMi - r,)r, (15) 



and 



v+n^p+e'^ 



POO 

AJ^ dp,plpMl-nf^. (16) 



Starting from A^_^p^^_|_^, and taking into account that Pedpe = E^dEf., Ey = E^ — Am and Ey = py, it is 



e 

possible to show that, using changes of variable, 

POO 



e +p— 



POO 

^ = A dpyplpeE,{l - r^)fl (17) 

•J 

In the regime x = E/T ^ 1, can be approximated by 

n = e"^' + ^^x^e-^' with Xi = E^/T. (18) 
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Since during the period of freezing, the temperature T is low compared with the energies appearing in 
the reaction rates, it is a good approximation to neglect Pauli factors. Then, to analyze the relationship 
between 

K-+p^n+. = ^ \ dvAeEJl (19) 



and 

\\ ^ =A 

^ >-i-'n -~^r\-A-(=> J 



roo 

A / dpAEJ^, (20) 
Jo 



it is enough to see how and relate themselves. From Eq. (Il8|), and recalling that in the reactions 
1/ + n ^ p + e~, El, = E(, — Am, we can directly write in terms of . This yields 



n = e-^-l^r, + e-^-/^^(2i?.Am + Am^)^^. (21) 



Then, the PDB is no longer valid, since the reaction rates in this case are related by 

K-^,^n^. = --^"^'^K^n^,^.- + ^e-^-/^ A dvAeE.ilE^^m + Am^)e-^^/^ (22) 

with Ey = E(. — Am. We shall say that we are working within a detailed balance with non-extensive 
corrections (BNE). 

It is worth noting that we can obtain an approximate expression among the rates, similar to what one 
obtains for the standard situation. Indeed, using Ey/T > 1 in Eq. (||), we get 

= e-^'"/^/,^ (23) 

what immediately yields 

AV ^ = e-^™/^AV ^ ■ (24) 

We can recover the standard PDB, invoking this additional approximation, what we shall refer as to stan- 
dard PDB (BST). In what follows, we shall show results both in the full BNE as in the BST approximation. 



and we shall discuss the range of validity of Eq. (23). 



3.2 e"^ + ri ^ p + z/ 
A similar analysis yields 

1 _ A 



^™/^A'p+P^n+e+ + V""^'"^''^ [ dp^plp.Ey{2KAm + Am^)e-^=/^, (25) 
where E^ = E^ + Am and we neglected the electron mass with respect to the energies. In fact, a careful 



look at Eqs. (22) and (Pq) show that both integrals are identical. 



Within BST, the relationship between the rates again simplifies, to give 

K++n^p+y = ^ ^"^^'^K+p^n+e+- (2^) 

3.3 n — > JO + e~ + 

Because of the fact that we have taken 1 — /* as 1 in Eq. (pT]), no modification arises in this case: 

- = Xn^p+e-+P = 0.0157^Am^ (27) 

T 
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This allow us to determine A in terms of a measurable quantity, the neutron mean life r, 



A = — Am^, 
r 4 



a = 255. 



(28) 



We shall neglect the free decay of the neutron in this section, particularly when computing A„j,. This 
makes Eq. (|8|) to assume the form A^^ = A^^^_^p_|_g_ + A^+^^p^-. Neglecting the blocking factors in Eqs. 
@ and (|lOp it is evident that this two terms are identical. The same happens with Eqs. (^) and (|13|), 
but in this case we have to neglect, in addition, the electron mass in order to have po equal to Am. Due 
to this, A9 = 2A'^^ ^ _ and A^ = 2A''__^ 

Using all previous results, the valid relationship between A^p and A^^, both in the BNE as in the BST 
approximation are 

g-Am/T^g^ + (g - l)e-^™/^/, (29) 



A'' 



where we have defined 



A 1-°^ 

with Ee = Ey + Am, for BNE, and 



/•CO 

Jo 

-Am/Tyi 



^np 



(30) 
(31) 



for the BST. 



I'+n— »p+e 



4 Explicit form for the rates 

To compute A^^ and A|„, it is clear that we have to know just one reaction rate, for instance A^, 
and the value of / defined in Eq. (]30|) . We can solve for A^^^^^^^. using Eq. (|20| ) . With t 
functions being as in Paper I we obtain 

roo 

A / dE^El{E^ + Amfr„. (32) 
Jo 



u+n-^p+e 



These integrals are easily computed changing variables to x = Ey/T and using the integral form of the 
function T{n). Noting that r(n + 1) = n! for n € Z, defining y = Am/T, and using (^), we obtain 



v+n^p+e 



1 a 
2ry5 



12 + 62/ + 



+ —{q-l) 



180 + my + 6y' 



These same considerations are valuable to compute /. It is given as 



i{y) 



1 a 

2ry5 



120?/ + 602/2 + 127/3 + 2/^ 



(33) 



(34) 



4.1 Xl^ 

We have already mentioned that A^„ = 2A 



u+n^p+e~ ' 
,2 



I.e. 



Kp{y) = ^ fl2 + 62/ + 2/'l + -^{q - 1) [l80 + 6O2/ + 62/^ 



(35) 



This equation is valid in both, the full BNE and the BST approximation. 
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Figure 1: Left: Standard reaction rates Ap„(y) and X^piv)- Right: Corrections produced by non-extensive 
statistics, normalized by (g — 1), using the full BNE. 



4.2 A"? 

We use Eqs. (|2|) and (|3|) to get, within the full BNE, 
It is also immediate to obtain this rate in the EST, i.e. 



120y + 60y2 + Uy^ + 



(36) 



(37) 



This last equation was used by Torres and Vucetich |p, where the BST was assumed from the beginning. 

It is worth noticing that since Ap^(y) and X^piv) have an asymptotic behavior analogous to that of the 
standard rates, see Fig. |l], also in this case we expect that in the limit of low temperatures (i.e. y oo) the 
neutron abundance behaves as a constant. That is, for low temperatures, the reaction rates are practically 
zero and then, processes interchanging neutrons into protons and viceversa are negligible. In this sense, 
X'^{y oo) stabilizes towards a constant. 



5 Evolution of the neutron abundance 

The formal solution for the evolution of the neutron abundance is given by Eq. (0). We shall use the 
variable y = Am/T, so that 
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and the factor / becomes 

=exp ("/J^^ (j) A^(y)^. (39) 

To evaluate the Jacobian dt/dy, we recall that the scale factor, R, in a Friedmann- Robertson- Walker 
metric, goes as i? ~ i/T, independently of the statistics Then, T/T = —R/R, and the rhs is given by 
Einstein equations, 

Here, pg is the energy density of relativistic species, given in Paper I. When the universe is dominated by 
e~, e"*", u and 7's, and then gh = 2, gj = 2 + 2 + 2 x ^ = 10 and g = J2b9^ + | = 43/4, we obtain 



We can now compute dt/dy, 

Since, t/T = -R/R, 

and to first order in (q — 1), 

In this sense, we see that 
where b and c are given 



Pg = ^gT' + 35.85{q-l)T\ (41) 
dt dt dT 1 Am 



(42) 



T (^-^Pg ) , (43) 



dy dT dy T y"^ 



1 _ 1 / 45 



15 

1 ^35.85(g - 1 



(44) 



^ = -&y[l-c(g-l)], (45) 
dy a 

/ 45 \ a 15 
6= -^7^ — r^' c= ^35.85. 46 

Numerical values for these constants can be obtained taking into account that: a = 255, g = 43/4 



Am = 1.29MeV, r = 887 ± 2s @ and G = m~? and are b = 0.25 and c = 5.07. This follows the treatment 



given in Ref. ||3|. 

Using all previous remarks, 

P{y,y')=e^v{K^{y)-K'^{y')), (47) 

with K'^{y) given by 

K'^{y) = -J dy(^]A'^{y), (48) 



or, equivalently, using Eq. (pSD , 

Tb 



K'{y) = --(1 - c{q - 1)) / dyyA'iy)- (49) 
a J 

We then see that the needed steps to solve for the evolution of X"^ as a function of y are: 1) find 
the function K'^{y) as in Eq. (|9|), 2) substitute this result in the expression (|^) for I'^{y,y'), and 3) 
finally solve the integral in (p^). The complexity of all functions involved requires a numerical procedure. 
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However, we can do further before going to numerical work, since our aim is to primarily see how the 
asymptotic value of X, X'^{y oo), is affected as a function of the parameter {q — 1), and ultimately, 
how this asymptotic value changes if one assume the full BNE as compared with the approximate BST 
approach. 



6 Standard Balance: BST 

In this case, X'^ {y) / M {y) coincides with its standard analogous X^^ (y) / A^^ (y) . Indeed, by definition, A'^(y) 
^pniy) + ^npiy)^ ^-iid since BST establishes A|„(?/) = e~^A^p(y), we obtain X'^ {y) / M {y) = (1 + e' 
X'\y)/K'K 

We can get K'^{y) taking into account that, using Eqs. (35) and p^). 



Computing the integral in Eq. (|49|), 



(l + e^^) 12 + 67/ + 2/2 +(g-l) 180 + 60?/ + 



=6(l-c(g-l))[Ki(y) + 



q-1 



K2iy)]. 



with 



and 



Kiiy) 



K2{y) = 

and where the functions Ei{l,y) are 



15 














+-) 


y2 


yJ 




yJ 



Eiihy) 



oo p-yt 



-dt. 



(50) 

(51) 
(52) 
(53) 

(54) 



We recognize in Eq. (51) the standard term K^^(y) = hKi{y). The presence of the factor 6(1 — c{q — 1)) in 
the Jacobian dt/dy produces the appearance in K2{y) of a second order term in (g — 1). It is important to 
note that because this function diverges as at the origin, it is not a priori obvious that we can retain 
(as was done in |Q]), only the term linear in {q — 1). This is a key remark, and we shall come back on this 
later. 

We shall then try to obtain the whole solution, without linearizing any of the functions involving 
{q — 1). Having shown that X'^ (y) / A'^ (y) = (1 + e^)~^, we can write, using Eq. (|^), the asymptotic value 
X'}{oo) = limy^^Xi{y), as 



X«(oo) = lim e^'(^) r 
y-"^ Jo 



jT^dy'- 



(55) 



This expression allows for an analytical study. Based on the divergent behavior of K'^ at the origin and 
because of the good behavior of 



d ( xuy'y 

dy' [ Ai{y') 



(1 + ev' 



(56) 



near zero, it is possible to discard a priori, i.e. without having the need to compare with any observation, 
the range of {q — l)-values such that the non-extensive theory, considered within the BST, looses its sense 
(i.e. the cases in which the asymptotic values predicted for X'^ are out of the interval [0,1/2]). Indeed, in 
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order to avoid problems at the origin in Eq. (^), K'^{y) has to be positive when y ^ 0. Using Eqs. (^T]), 
(|5^ ) and (|5^), this condition translates immediately into an allowed interval for {q — 1) 



smce 



lim K'i{y) ~ lim 6(1 - c{q - 1)) 



y3 ^ y3^ 



+ 



30 30 „ 

^ + — e^J' 

y 



I.e. 



lim K«(?/) ~ b{l - c{q - 1)) 



+ ('7-1 



30 



(57) 



(58) 



_r r 

We see that, to have K'^{y) > 0, (g — 1) must be such that the inequality —8/30 < (g — 1) < 1/c is fulfilled. 
Recalling that c = 5.07 we see that values that do not fulfill the condition 



0.27 <{q-l) < 0.2 



(59) 



must be automatically discarded. We mention in addition that no matter the sign of K'^, the factor e^^^^ 
multiplying the integral in Eq. (^) tends to 1 when y ^ oo. 

Fig. shows the behavior of the integrand of Eq. (p^) for different (q — l)-values. Note the abrupt 



change of the plots when (q — 1) is near the extremes of the range given by Eq. (59). The analytical 
conclusion is then reinforced by these plots and is completely confirmed when numerical computations are 



made. These latter show that the integral (55) grows without limit when (q—l) assumes values out of the 
range given by Eq. ([5^). We have then showed that there exist a range of {q — 1) such that the asymptotic 
value of X'^(oo) has physical sense. We can now ask if there also exist a range of {q — l) G [—0.27, 0.2] such 
that the X'^{oo) obtained to first order in (q — l) is consistent with the whole computation. This would be 
so if we can prove that the first order result do not differ much from the real result that we have already 
got. 

6.1 First order computation 

To compute X'^{oo) to first order in (q — l), we need to write Eq. ( p5[ ) to first order and neglect quadratic 
terms in Eq. (|5ll). This gives 

K^{y) = K'\y) + {q-l)K'{y), 

with 

K'^\y) = b 



and 



We write e^'' to first order as 



4 3 1\ / 4 1 > 

— + — + - + — + — )e y 

yi yZ yJ \^yi y2 



30 15 3\ /30 3 , „ ^ , , 
— + — + - + - + - e-^- 3Et{l,y) 

r y yJ y 



e^' = e^'e^-^-i)^^ ^ e^"[l + {q - 1)K% 



cK 



St 



Substituting this expansion into (^5|) and retaining only the linear term we get 

X9(oo) = X-'\oo) + {q- l)X^(oo), 

with 



and 



X-'\oo)= lim e^"(j^) f\-K'Hy')^^dy', 
^ ' Jo {l + ev')^ 

X^(oo)= lim e^'*(^) [\-^^'^y'\K\y) - K\y')]—^^dy' . 



(60) 

(61) 

(62) 

(63) 

(64) 
(65) 
(66) 
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y y y 

Figure 2: Behavior of the integrand needed in the computation of X*(oo) for different values of (g 
within the BST approach. 
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Figure 3: Comparison between X'^(oo) as a function of — 1), to first order, and using the full numerical 
solution, both obtained within the BST approach. 



This integrals can be computed numerically, and the result is 

X''(oo) = 0.15 + (g- 1)0.18 (67) 
where the value X'^*(oo) = 0.15 is the standard one. Fig. |^ shows the linear dependence of Eq. (|67|), 



together with the full solution given by Eq. (55). Although within the framework of the asymptotic 
approach to quantum distribution functions, an analogous equation to ( |67| ) was presented already by some 
of us 1^], we have now shown, for the first time, what is behind that approach, what sustains it, and the 
range of its validity. 



7 Non-extensive balance: BNE 

In this section we shall solve Eq. ( [3^ ) without approximations. That is, we shall take for and A^p, the 
expressions given by ( ^5| ) and (|36|), respectively. Our aim will be to quantify the amount of the deviation 
that it is produced when the BST, instead of the full BNE, is considered. 

By definition, A'^{y) = Xf,n{y) + A^p(y), adding (|35|) and (|36|) we see that A'?(y) is given by the sum of 
two terms 



A^*(2/) 



-(12 + 6y + y2) [l + e"?'] , 



1 a 



A^{y) = 90 [1 + e-y] + 30y [1 + Be-v] + y^l + 21e-y] + 12y^ + 



(68) 
(69) 



Because of the linearity of the integral operator, the same happens to K'^{y) in Eq. (p9|). Using the 
previous expressions, we get 

1 



K'i{y)=h{l-c{q-l))[K^{y) + 



-K2{y)]. 



(70) 
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with 



and 



Ki{y) 



4 3 1 \ / 4 1 , „ 

— + — + - + — + — )e y 

y.i yZ yJ ^ yd yZ 



K2{y) 



30 15 3\ /30 60 3 , „ , 

y6 yZ yJ ^yd yZ y I 



(71) 



(72) 



where Ei{l,y) were given in (54). As before, we do recognize the standard term K^^(y) = bKi{y). Also 
here, within the BNE, it happens that 



hm 



0. 



(73) 



Indeed, since within the BNE, in the hmit y ^ oo we have A^„(y) oc y~^ and X^piv) °^ y~^^~^ ■> we obtain 



hm 



hm 



hm 



1 



yix^ ^q^y^ jP^ci) 1 + A^n(y)/A^(y) y^'^ l + 
We then obtain an analogous to Eq. (EEi) 



0. 



(74) 



X^(oo) = - lim e^'^-y^ [\ 
y^<^ Jo 



-K'>{y')J_ ( ^npjyY 

dy' \ \ln{y) 



dy'. 



(75) 



It is important to note that K'^{y) diverges as y~^ at the origin, and as —y at infinity. Then, first order 
developments will not do. The complex dependencies of the integrand of Eq. ( |38t) with {q — 1) makes much 
harder to a priori analyze the validity range, as done within the EST. We can, however, make a detailed 
analysis of the behavior of the derivative of the function \^„{y)/M{y) and of the integrand of Eq. (^8|), 
for different {q — l)-values. We show this in Figs. ^ and g. Prom the analysis of these figures we see 
that the EST and the BNE differ in a fundamental way, which translates into the value of X'^(oo). Within 
the BNE, and because of the behavior of the derivative of the function X^^{y) / A'^ (y) for negative values 
of (g — 1), it is not guaranteed that the integrand in Eq. ( ^8] ) is negative. For (q — 1) < —0.1 the integral 
of ( ^8|) is positive (i.e. X'^{oo) < 0), and this discards values of q such that {q — 1) < —0.1. On the other 
hand, looking at Fig. |^ for positive values of (q — l), we see that the upper bound is between 0.19 and 0.2, 
in this case similarly to the BST. Then, for the BNE, the values for which the theory has physical sense a 
priori of any experimental comparison are 

- 0.1 <{q-l) < 0.2. (76) 
The BNE is more restrictive than the BST. Out of this range, the physical sense of the description is lost. 



8 Comparing both approaches, and bounding the value of q 

We now compare the results obtained in Sections ^ and Fig. ^ shows X'^{oo) for different values of 
{q — l), both in the full BNE as well as in the BST approximation. Differences are notable, and they grow 
with increasing values of Ig — 1|, as it tends to its limiting values. Table 1 gives the explicit values of A'^(oo) 
obtained using a numerical integration scheme for the BST and for the full BNE, and the analytical results 
of the first order approximation for the BST. For values of {q — 1) near 0, the BST numerical computation 
is perfectly described by the linear approximation. The deviation from the standard case is greater for the 
BNE than for the BST, for equal values of {q — l), i.e. the BST underestimates the effect of non-extensivity. 
This can be seen directly from Fig. ^, or equivalently from Table 2, where we show the coefficients of the 
linear fit of X'^{oo) near 0, both for the full BNE and the BST. In the sake of completeness, we also show 
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(q-1)=-Q.2 



(q-1)=-Q.17 



1Q 15 



1Q 15 





(q-1)=-Q.1 



1Q 15 2Q 



standard 



1Q 15 2Q 




(q-1)=0.1 



1Q 15 2Q 



(q-1)=0.15 



1Q 15 



(q-1)=0.19 



1Q 15 



(q-1)=0.2 



1Q 15 



;ure 4: Behavior of the derivative of X^^{y)/A^{y) for different (g — l)-values within the BNE. 
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Table 1: Comparison of X*(cxD)-values as a function of (g — 1) for the different approaches. 





XI (BST) 


X'i (BST) 
(0 (g-1)) 


X'i (BNE) 


-0.150 


0.1472 


0.1230 


-0.6571 


-0.120 


0.1441 


0.1284 


-0.1061 


-0.100 


0.1431 


0.1320 


0.0261 


-0.050 


0.1440 


0.1410 


0.1299 


-0.010 


0.1484 


0.1482 


0.1476 


-0.001 


0.1499 


0.1498 


0.1498 


0.0 


0.1500 


0.1500 


0.1500 


0.001 


0.1502 


0.1502 


0.1503 


0.010 


0.1520 


0.1518 


0.1523 


0.050 


0.1625 


0.1590 


0.1635 


0.100 


0.1851 


0.1680 


0.1915 


0.120 


0.1991 


0.1716 


0.2109 


0.150 


0.2303 


0.1770 


0.2561 
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Table 2: Coefficients of the linear fits of X'^{oo) near (q — l) = for the full BNE, the BST approximation, 
and for the first order computation of the BST. 



Case 






BST Approximation 
BST 0{q-l) 
BNE Full 


0.15 
0.15 
0.15 


0.18 
0.18 
0.23 



the ffi'st order result for the BST (this is not a fit, but comes directly from the analytical computations 
shown before). 

Until now, we have discussed the physical viability of the statistical description, establishing a range 
of a priori discarded values of {q — 1) both for the full BNE case, the real situation, and for the BST 
approximation. These values are those for which X'^{oo) is < or > 0.5. This can be done before than the 
comparison with any given observational data or experiment. Recall that X(oo) is the neutron to proton 
ratio, for which clearly, a negative value has no sense. In addition, X{oo) > 0.5 is in conflict with the 
reaction rates dependence with temperature, since X = 0.5 is the initial condition for the kinetic equation 

when T oo (i.e. 7/^0). Within the range of physical viability we now look for consistency with 
observations. This will further restrict the range of admitted values of (g — 1). To obtain a direct bound 
upon (q — 1) using the primordial abundance of ^He it is necessary to study in detail the free decay of 
neutrons, happening between the moment of freezing out of the weak interactions (t ~ 1 s) and the moment 
in which the temperature of the universe is similar to the binding energy of D (t ~ 3 minutes). We shall 
provide, for the first time in a non-extensive setting, a detailed account of this in the following section, but 
nevertheless, let us give here some preliminary considerations in the sense of Ref. . 

Including the effects of the neutron decay in the equation for the evolution of X'^{t) we have, 

X'^{t)=exp{-t/T)X'^{t), (77) 

where X"^ is the already obtained ratio and r is the neutron mean life. In the capture time, t = tc, when 
the temperature is below the D binding energy (2.23 MeV), neutrons are captured into D. Then, almost 
all neutrons present at t = tc are converted into ^He. Substituting the value of tc into Eq. (|77|) and using 
the value of X'^{oo) previously found, we shall obtain half of the mass fraction of all ^He primordially 
produced. For the time being, we shall adopt the standard value for exp(tc/T), which is ~ 0.8 |^]. With 
this we obtain, 

= 2X1 ~ e-*^/^2X«(oo) 

= 0.8x2x [0.15 + (g- 1)0.23] 

= 0.24 + (g - 1)0.37. (78) 
There is no absolute consensus about the observational value of '^He. The two greatest compilations give 



m 



Yp = 0.244 ± 0.004 and Yp = 0.234 ± 0.004. (79) 
To consider a typical case, we average over the two mean values and twice the error bar, 

Yp = 0.239 ± 0.008. (80) 

If at the same time we neglect the difference between the standard theoretical and the observational 



values, which is given by 0.001, in order to obtain a first bound, we can get, using Eq. (78) and asking for 
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|g- 110.37 < 0.008, 

|g-l|< 0.021. (81) 

Within the most complete treatment of the principle of detailed balance (BNE), which accounts for the 
more in depth analytical study on the influence of a slight non-extensivity in primordial nucleosynthesis, 
the value of q could differ from 1 at the level of a few percent and still be in agreement with current 
constraints. See below for a more detailed treatment. 



9 The capture time 

The aim of this section is to show how the capture time tc is modified in the non-extensive framework 
Recalling that at early times R/R = {St: pq / Snipi)'^ and R{t)Ty{t) = constant after the e^e~ annihilation, 
it is straightforward to write 

3 \V2 /-oo 



t 



8ttG 



r + to, 



(82) 



where to is a constant whose standard value is of the order 2 seconds, see In Paper I, it was shown that 
the energy density could be written as pR = {-k"^ /30)g'iT'^. Taking also into account that the photon and 
neutrino temperatures are related by = (11/4)^/3 [i + 0.109(g - 1)]T^, we can solve the definite integral 
to obtain: 

t=(^ 



1/2 ,^ 2/3 



^) ^[1- 5.06(^-1)]+ to, 



4 J 



rp2 



(83) 



where Qeff = 2(11/4)^^/3 + 21/4 ~ 12.95. 

In equilibrium, neutrons, protons and deuterons behave as free non-relativistic gases, with number 
densities given by 



niiT^\^/'^ 
2-K 



-{mi-^ii)/T^i 



1 + 



q-l / 15 



+ 3 



nii - Hi 



+ 



mi - Pi 



(84) 



Here, the subscript i stands for n,p and D, for neutrons, protons and deuterons respectively, gn = dp = 
2,gD = 3 and Pi,nii are the chemical potentials and masses of the particles. At early times, these gases 
are in chemical equilibrium, so that pd = Pn + Pp- Therefore, 



'n'qn^ _ 9n9p firinmp 



n- 



D 



9D 



nriD 



3/2 /TH X 3/2 
22. \ 
2tt 



u{{mn - Pn)/T)u{{nip - Pp)/T) 



u{{niD - Pd)/T) 



(85) 



with eo = + nap — tud, and where we have made use of the definition. 



u{{m^ - pi)/T) = 1 + 



1 / 15 



rrii - Pi 
T 



+ 



mj - Pi 
T 



(86) 



Equation ( |85D is the generalized Saha equation that describes the deuterium formation (see Paper I for 
a general derivation). Helium appearance is inhibited by the "Deuteron bottleneck", represented by the 
reactions: 



n + p 
D + D 
D + T 



D + 7, 
T + p, 
te + n. 



(87) 
(88) 
(89) 
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Eq. (p^) can be given in terms of the normalized abundance fractions X* = n^^/n^ . To this purpose we 
introduce the baryon to photon ratio r/g = /n^ and recah that = [2^(3) + 3!C(4)(g — l)](T^/7r^). We 



get for the quantity = X^X^ /Xj^ , to first order in {q — 1), the following result 



G1 ^ ^ J_ fZV!^y^%-.z,/T- 



q — 1 , s 
1 H ^(^^ + Up- ud) 



(90) 



Here, as in Paper I, Ui = 15/4 + 3(mj — Hi)/T + ((mj — ^i)/T) . As in the standard case analyzed in Ref. 
1^], to determine the neutron capture time tc we need to examine the sequence of reactions (87-^). It is 
convenient to do so in terms of an scaled temperature z = eo/T^f, so that they involve rate parameters of 
the form 

i?'^ = #MTnf, (91) 



where {crv)^ denotes the thermal average of the relevant cross section times the relative velocity Using 
Eq. (83) and writing the baryon number density in terms of the baryon to photon number density ratio, 
r]g, we get 

= I (^^) i^T' C{5)eDm,i{av)T[l + 0.34(g - 1)]. (92) 
Taking r]Q = 3.57 x 10^'^ as a nominal value for rjg and considering {anpv)T = 4.6 x lO"^'^ [^] we obtain 

Rl^ = 3.63 (^^y ^[1 + 0.34(g - 1)]. (93) 

Neutron and proton populations are mainly determined by the kinetic equations 
dX^ dX^ 

— —R^^plXgX^ — C^^pXg ], — ~^np[-^q-^q ~ ^np-^q ]• (94) 



As in the usual case, if the neutron population is not depleted by other reactions such as Eqs. (|88- 
89), protons, neutrons and deuterons are kept in equilibrium, relationships X^ + X'^ + 2X'^ = 1 and 
X'j^ = (G^p)~^X^X^ being valid. Since (G^p)"^ is very small, the deuterium number density will be also 
small, and we can write a first approximation as 

^.""'^'^ = {Glp)-^X-/'^Xl'^'\ (95) 

where Xg'^^^ and Xg'^^'^ are the unperturbed populations, that obey Xg'^^'^ + Xg'^^^ = 1. Using this we 
now have 

X^p+X?,c,l- 2{Glp)-'X^'('^XPg'('\ (96) 

Recalling Eq. (|90|), one sees that for the z- values we are interested {z ~ 30), the major z dependence of 
{G^p)~^ goes as e^. Hence, to first order, 

^^{X^p + X^J ^ -2(G^,)-iX,"'Wx,^'W. (97) 



Adding Eqs. (|94|) and using Eq. (|95D we find that 

= Rlp[X^XP - GlpX^]. (98) 
When deuterium number density is depleted, its number density will change according to 



dX 



D 



= Rlp[X-gXl - GlpX^] - RUm^f - GUxJX^p] + ... (99) 
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where Rdd the scaled rate for the reaction ( |88| ) and G^j^, is the Saha factor which gives the equiUbrium 
value for the ratio {X'^)'^ / Xj . It can be shown that, in contrast to Gnp, Gdd is always a small number 
(see Section III of Ref. Q). In fact, in the standard framework Gdd is of order e~^^. We shall obviously 
be safe in considering G'j-,^ = Gdd, and therefore in neglecting it. We can obtain R%d inserting the 
value of {(Tddv)t in Eq. (|92[). The calculation of {c7ddv)t is a non-trivial task, which involves the use of 
a phenomenological fit to the cross section plus a thermal average carried out using Boltzmann velocity 
distributions for the deuterons. For the sake of simplicity, we shall adopt the value of the thermal averaged 
cross section used by Bernstein et al. This can be written as 

2/3 

{(Jddv)t = 0.016 ^ e-^-^^'"\ (100) 
Introducing the nominal value of rjo as we did it when we wrote Rnp, we can now write for Rdd 

R%^ = 1.76 X 10^ (^) ^-4/3g-1.44.V3 ^ _ ^^Q^^ 

When the deuterium number density decreases enough, the chain of reactions that converts almost all 
neutrons into helium is initiated. We may therefore identify the temperature T^,c, at which the neutrons 
are captured, or equivalently, the value of z = Zc = ezj/T^^c, as in the standard case, by the condition that 



'dxl^ 

dz 



0. (102) 



Z = Zc 



Once we have neglected the factor G^j^j, the condition stated above, together with Eqs. (p5| ) and ( |98[ ) and 
the approximation X^ = Xq'^'^\ gives 

2Xf ^ 1. (103) 
In order to solve this equation we proceed as follows. Let us first rewrite the previous equation as 

Glp = 2Xf')xl'^'^Rl^. (104) 

Note that we can put the rhs in terms of z. In Section 8 we have already given the asymptotic value of Xq'^^^ 
as 0.15 + 0.23(g — 1), and therefore the value of Xq'^^'^ can be found trough the relation Xq'^^'^ + Xq'^^^ = 1. 



Now, using Eq.( |10lD leads, to first order in (g — 1), to 

Glp = 4.5 X 10^ (^^^ [1 + (g - l)1.6]z-4/3g-i.44.V3_ ^^Q^^ 

On the other hand, we also have an expression for G^„: it is given by Eq. (TO). Therefore, if we were are 



able to write Eq. (|90| ) in terms of z we can equal both expressions and solve for the unknown Zc, which 
will satisfy the equality for a given value of {q — 1). Unfortunately, this is not easy. The difficulty is due 
to the presence of the chemical potentials in in combinations other than the relation /xd = + Aip- 
In this sense, we confront here a problem similar to that found in the recombination study of Paper I. We 
shall just mention that it is possible to show that the function Un + Up — ud has 10^ as an upper bound, 
see Appendix for details. Therefore, in the range of z-values that concerns us, the corrections in G'^p due 
to non-extensivity will be at most of order [{q — 1)/2]10^ and we can then write for G'^p = G'^p{z), 

Glp = 2.98 X 10^2 ( ^ ) [1 + 5 X W\q - l)]z^/^e-\ (106) 
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Figure 7: Left: as a function of (g — 1). Right: Capture time as a function of (g — 1). Both plots are 
constructed for a nominal value of riq = 3.57 x 10^^'^. Note the asymmetric range of the (g — l)-values 
considered. See text for explanation. 



Thus, equating expressions ( |105D and ( 106 ), and taking ijq = rjQ, we finally found that Zc must satisfy the 
following relation. 

The standard result (i.e. {q — 1) = 0) gives Zc = 27.08. The numerical solution of this equation for 
different values of the parameter {q — 1), and the corresponding capture time, are shown in Fig. It is 
interesting to notice that while values of (g — 1) > can be used within the analytical study presented in 
this Section without any limit, this is not the case for values of the parameter such that {q — 1) < 0. The 
problem arises in that in these cases, Gfip can become negative, what clearly has no sense. In particular, 
for (q — 1) = —0.002, the denominator in Eq. (107) is identically zero, while for values of (g — 1) < —0.002 
it is negative. Although this mines the physical viability of the description for these values of q, in the same 
sense that what was referred before, other assumptions are taken in this Section, and only a numerical 
code can give a more precise answer in either way. Note that to obtain Fig. 7 we have made use of the 
non-extensive time-temperature correction. In Fig. 8 (left panel), we show the dependence of the ^He 
abundance as a function of {q — 1), taking into the free neutron decay correction. 

Some comments about Figs. 7 and 8 are worth doing. Firstly, note that the dependence of the ^He 
abundance, when neutron free decay is taken into account, is no longer linear as it was before, see for 
instance Eq. (|7^. This makes the study of the nucleosynthesis process much more interesting, and the 
use of the full numerical code a worth doing task. Now we know that even when considering a first (linear) 
order correction in the quantum distribution functions, the final output in the primordial nucleosynthesis 
cannot be reduced to a linear correction in the abundances. Secondly, it appears that for a given value of 
Yp (and ij) there are two values of {q — 1) that could fit well the observations. Finally note that the y-axis 
of Fig. 8 (left panel) is showing variations at the level of less than 1%, and in principle, all values of {q — 1) 
shown there are admitted by the current constraints on the primordial abundance of ^He. 

In Fig. 8 (right panel) we show for the first time the dependence of the ''He abundance as a function 
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Figure 8: Left: ^He abundance as a function of {q — 1), taking into account neutron decay corrections. 
This plot is constructed for a nominal value of rjq = 3.57 x 10"^''. Right: Primordial abundance of ^He as 
a function of r]q. Several {q — l)-values are shown. Two horizontal lines stand for a conservative current 
constraint on Y„. 



of r]q, for different {q — l)-values. The general trend of an increasing abundance as a function of r] is 
maintained, but for different values of q, the curves cross the current constraints in different positions. 
This plot actually represents a 3D problem, with a two parameter space. The influence of t] appears when 
numerically solving the equation for Zc- Indeed, the computation of the correction to the value of Yp in the 
full analytical treatment presented here has to take into account different contributions. A direct linear 
dependence on (g — 1) is within the time-temperature relationship (see the first equations in this section, 
especially Eq. (|83|)) and in Eq. (|78|). In addition, we can see that if Zc is smaller than in the standard 
case, then Tc will be bigger, and tc will be smaller than their respective standard counterparts. Thus, e^*/"^ 
will in turn be bigger, what happens for instance in the case of a negative value of {q — 1), and this would 
provide a bigger Yp. However, the linear correction in Eq. ( |78| ) is competing with that positive deviation; 
this is what finally makes the curve to deviate from a straight line. 



10 Conclusions 

In this paper, based on previous works (Paper I), we have revisited the problem of the primordial genesis 
of light elements within the framework of non-extensive statistics. We have particularly paid attention to 
the form of the principle of detailed balance that is valid in this new setting. We have shown that its usual 
form is no longer valid, but instead, that it is just an approximation, dubbed here as the BST. The full, new 
form of the principle of detailed balance was derived and named as BNE, and we have studied it in detail 
too. By doing so, we were able to disguise the range of validity of previous approaches to the topic, and 
to formally see the origin and range of the deviations between the standard and non-extensive scenarios 
for the primordial production of elements. We have also analyzed the neutron free decay correction to 
the capture time, and found that even when considering a first (linear) order correction in the quantum 
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distribution functions, the final output cannot be reduced to a linear correction in the abundances, as Fig. 
8 explicitly show. By comparing with the latest observational data, we have obtained a new bound on the 
upper limit of the parameter |g — 1|. We think that the bound being presented here (referring also to Figs. 
7 and 8) is obtained within a much more detailed basis than the previous ones, see for instance [||, 
and then that this bound should be considered as the more reliable of them all. Indeed, we believe that 
the only way to improve it is to go to a direct implementation of the primordial nucleosynthesis code. By 
doing so we would have the possibility of studying the modifications to the abundances of all light elements 
in a simultaneous way. This further extension is currently under analysis. 

Appendix: Un + Up — ud 

To have a rough idea of the order of magnitude of the function Un+Up — U£), we need to compute {rrii — ^i) /T 
in the range of temperatures we are interested in, with i = n,p, and D. Firstly, we note that for this purpose 
it is enough to consider the standard values of all these quantities. This is because non-extensive effects 
here would introduce only second order corrections in (q — 1), that we are disregarding in our computation. 
Recalling the relation Xi = {gi/nB){miT /2tt)'^^'^ exp{{mi — fii)/T), we see that in order to estimate the 
values of (mj — fJ'i)/T we need to know the abundances Xi. Then, to get Un + Up — u^, we have made the 
following steps: 

1. Using the time-temperature relation we can put Xn{t) = exp{—T/t)Xn{tfrezee) ^ an explicit function 
of T. Thus Xn{T) follows, as well as (m„ — fin)/T. 

2. Considering Xn{T) together with the obtainable function X^Xp/X^ (remember that fJ'n + fJ-p — fJ-D = 
0), and imposing the (approximate) constrain Xn + Xp + 2XD = 1, we get the values of Xp and X^. 
Then, the corresponding values of {nip — Hp)/T and (m^ — ij,d)/T can be obtained. 

These results allow us to construct Un + Up — ud and thus to see that 10^ constitutes an upper bound in 
the range O.OTMeV < T < O.lMeV (i.e. 22 < z < 32). 
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